Ineffectiveness of Pade resummation techniques in post-Newtonian approximations 
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We test the resummation techniques used in developing Pade and Effective One Body (EOB) 
waveforms for gravitational wave detection. Convergence tests show that Pade approximants of the 
gravitational wave energy flux do not accelerate the convergence of the standard Taylor approxi- 
mants even in the test mass limit, and there is no reason why Pade transformations should help in 
estimating parameters better in data analysis. Moreover, adding a pole to the flux seems unnec- 
essary in the construction of these Pade-approximated flux formulas. Pade approximants may be 
useful in suggesting the form of fltting formulas. We compare a 15-orbit numerical waveform of the 
Caltech-Cornell group to the suggested Pade waveforms of Damour et al. in the equal mass, non- 
spinning quasi-circular case. The comparison suggests that the Pade waveforms do not agree better 
with the numerical waveform than the standard Taylor based waveforms. Based on this result, we 
design a simple EOB model by modifiying the ET EOB model of Buonanno et al., using the Taylor 
series of the flux with an unknown parameter at the fourth post-Newtonian order that we fit for. 
This simple EOB model generates a waveform having a phase difference of only 0.002 radians with 
the numerical waveform, much smaller than 0.04 radians the phase uncertainty in the numerical 
data itself. An EOB Hamiltonian can make use of a Pade transformation in its construction, but 
this is the only place Pade transformations seem useful. 

PACS numbers: 04.25.dg, 04.25.Nx, 04.30.Db 
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I. INTRODUCTION 

Even though general relativity was developed at the 
beginning of the twentieth century, no analytical solution 
is known for the two-body problem. Until recently, at- 
tempts to find a numerical solution failed because of the 
complexity of the mathematical equations and the insta- 
bilities inherent in the analytical formulations being used. 
In thepast few years, breakthroughs in numerical rela- 
tivity [J 13, [1, 01 allowed a system of two inspiraling black 
holes to be evolved through merger and the ringdown of 
the remnant black hole [5, 6, 7, 8, 9, 10, 11, 12, 13]. 

Studying the late dynamical evolution of these inspi- 
raling compact binaries is important because they are 
among the most promising source of gravitational waves 
for the network of laser interferometric detectors such 
as LIGO and VIRGO. The detection of these gravita- 
tional waveforms (GW) is important for testing general 
relativity in the strong field limit. Moreover, these de- 
tectors can extract from the waves physical data about 
these sources such as the component masses and spins 
and the orbital eccentricity. For an unbiased extraction 
of these parameters, a large bank of accurate waveforms 
needs to be constructed. Numerical relativity alone can- 
not compute all the waveforms needed because of the 
computational cost. Instead, the waveforms are based 
on post- Newtonian (PN) approximations [3, [lB| ■ 

The post-Newtonian approximation is a slow-motion, 
weak-field approximation to general relativity. In order 
to produce a post-Newtonian waveform, the PN equa- 
tions of motion of the binary are solved to yield explicit 
expressions for the accelerations of each body in terms 
of the binary's orbital frequency [H, [H, |l3, IH, IH, 
H HI m, ii, m ill . Then solving the post-Newtonian 



wave generation problem yields expressions for the gravi- 
tational waveform h and the gravitational wave flux F in 
terms of radiative multipole moments 26] . These radia- 
tive multipole moments are in turn related to the source 
multipole moments, which can be given in terms of the 
relative position and relative velocity of the binary [2^ . 
Instead of comparing the post-Newtonian waveform with 
a numerical waveform along a certain direction with re- 
spect to the source, the comparison can be done in all di- 
rections by decomposing the waveform in terms of spher- 
ical harmonic modes. For an equal-mass non-spinning 
binary, the (2, 2) mode /122 [H, [13, [Mi is often used 
to compare numerical and post-Newtonian waveforms be- 
cause it is the dominant mode. Its time derivative /122 is 
used to compute the gravitational wave flux. The result- 
ing expressions for the orbital energy E, the gravitational 
energy flux F and the amplitude /i22 are given as Taylor 
series of the frequency-related parameter 
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where M is the total mass of the binary and G = c = 1. 
The invariantly deflned "velocity" 



(2) 



another dimensionless parameter, is often used in writing 
these Taylor series. 

Computing PN series to high order is difficult and time 
consuming. Since the various PN expressions are given 
as sl owly convergent Taylor series, the Pade transforma- 
tion [22, [s^l was suggested in Ref. [13] to accelerate the 
convergence of these series. The Pade transformation, 
P™, consists of writing a Taylor series, Tk, of order k as 
the ratio of two polynomials, one of order m in the nu- 
merator, and another of order n in the denominator, such 
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1.0000000 


1.0000000 
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3.0099999 


3.0099999 
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5.0300499 


-401.0000 


3 


6.3834834 


9.1313636 


4 


7.0635838 


7.0601492 


5 


7.3369841 


7.4053299 


6 


7.4285732 


7.4747817 


7 


7.4548724 


7.4645660 


8 


7.4614801 


7.4631404 


9 


7.4629558 


7.4633014 


10 


7.4632524 


7.4633191 


11 


7.4633066 


7.4633174 



TABLE I: Convergence of the Taylor expansion, Exp^ — 
X]fe=o 'f'' /k\ of the exponential function Exp(u) and its Pade 
approximant Exp™^'' at u = 2.01, m = [n/2j. The Pade ap- 
proximant converges to six significant figures while the Taylor 
series converges to four significant figures ai v — 2.01. The 
error between the exact value of the exponential, 7.46331734, 
and the Pade approximant Expg(ii = 2.01) is 6 x 10~*, while 
the error between the Taylor approximant Exp^^{v = 2.01) 
and the exact value is 10~^. 
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FIG. 1: Convergence of the Taylor expansion, Exp„ = 
X]fc=o '^''Z'^' exponential function Exp(w) and its Pade 

approximant Exp^^^" at v ^ 2.01, m = [n/2]. The Pade 
approximant converges faster than the Taylor series. 



that m + n < k. If well-behaved, this method accelerates 
the convergence of a Taylor series as the order of the Pade 
transformation, m + n, is increased. For example, in Ta- 
ble U we compare the convergence of the Taylor expansion 
of the exponential function, Exp„(w)(= e"), at order n 
to its Pade approximant, Exp^'^''{v) = P™+''[Exp„(z))], 
along the diagonal, where m = [n/2j and e = or 1. 
After twelve terms {n = 11), the last two partial sums 
of the Taylor expansion converge to 4 significant figures. 
However, the last two Pade approximants Exp5(i;) and 
Exp5(w) converge to 6 significant figures. The error be- 
tween the exact value of the exponential, 7.46331734, and 
the Pade approximant Exp5(i; = 2.01) is 6 x 10~^, while 
the error between the eleventh order partial sum and the 
exact value is 10~^. Fig. [1] shows the convergence of the 
Taylor expansion of the exponential function and its Pade 
approximant. 

The hope of accelerating the convergence of the post- 
Newtonian Taylor series of the energy and flux motivated 
the use of their Pade approximants to construct Pade 
approximant waveforms O, [H, [H, [H H [H [13, SH, 
Ej S, IHj Ell ■ If these resummation techniques ac- 
celerate the convergence of the Taylor series in PN ap- 
proximations, the range of validity of PN approximations 
suggested by Ref 147] could be extended. Moreover, the 
work of Refs. [48|, \^ in the test mass limit motivated 
the addition of a simple pole to the flux F a, binary 
system as the bodies approach the light ring orbit. By 
mathematical continuity, the existence of a pole in the 
equal mass case was anticipated 34] . 

More recently, waveforms are constructed by including 
these ideas in Effective One Body (FOB) models. The 
FOB approach [H, HE [H, [H, M, M, ES, S, El M, 



IsU . m, m, m, m, [131 aims at providing an accurate 
analytical description of the motion and radiation of coa- 
lescing binary black holes. The approach consists of three 
separate ingredients: 1) a description of the conservative 
Hamiltonian part of the dynamics H, 2) a formulation 
of the radiation reaction force J- from the radiated flux 
F and 3) an expression of the GW waveform amplitude 
emitted by the coalescing binary system (i.e /i22)- 

The flux plays an important role in approxima ting th e 
radiation reaction force in the FOB models [4l','58','59|. 
The leading-order radiation reaction force [60, 61, 6^ 
enters the equations of motion at 2.5PN order. Since 
the equations of motion are known only to 3.5PN order, 
one has to rely on the assumed balance between energy 
loss in the system and radiated flux at infinity (63l . [64| 
to generate an approximate expression of the radiation 
reaction force at 3.5PN order beyond the leading term. 

Ref. [gHI computes the GW energy flux and GW fre- 
quency derivative from a highly accurate numerical sim- 
ulation of an equal-mass, non-spinning black hole bi- 
nary. By assuming energy balance, the (derivative of 
the) center-of-mass energy is estimated. These quan- 
tities are then compared to the numerical values using 
various Taylor, Pade and FOB models. The main goal 
of Ref. [63] is taking a set of well-established proposals 
in the literature for approximating waveforms and seeing 
how well they work in practice. Another goal of Ref. |63 
is to examine some modifications of those proposals. The 
main goal of this paper, by contrast, is to show that a 
key ingredient in those proposals does not appear to be 
necessary. 

In Ref. :66], Blanchet gave an argument that Pade 
and FOB resummations are unjustified because for two 
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comparable-mass bodies there is no equivalent of the 
Schwarzschild light-ring orbit at the radius r = 3M. His 
argument is based on the PN coefficients of the binary's 
energy and their relation to predicting the innermost cir- 
cular orbit. He finds that the radius of convergence of the 
PN series, which is related to the radius of the light-ring 
orbit, is around one (instead of 1/3 as for Schwarzschild). 
Blanchet concluded that Taylor series converge well for 
equal masses and that templates based on Pade/EOB are 
not justified because the dynamics of two bodies in Gen- 
eral Relativity does not appear as a small " deformation" 
of the motion of a test particle in Schwarzschild. This 
paper arrives at similar conclusions but not by consid- 
ering the innermost circular orbit, which is not precisely 
defined in the full nonlinear case. Instead, we compare 
Pade approximants of the flux and Pade/EOB waveforms 
to the numerical data of Refs. [btI. [6^. 

In this paper, we focus on testing two main techniques 
involved in building EOB models: the systematic use of 
Pade approximants, and the addition of a pole to the 
flux. The goal is to simplify these models by removing 
any unnecessary procedures in designing waveforms that 
provide good agreement with numerical waveforms. 

Damour et al. [H, HI] first suggested techniques for 
resumming the Taylor expansions of the energy and flux 
functions. Starting from the PN expansions of the en- 
ergy E and the flux F, they proposed a new class of 
waveforms called P-approximants, based on three essen- 
tial ingredients. The first step is the introduction of 
new energy-type (Eq. and flux type (Eq. I16p func- 
tions, called e{v) and /(u) respectively. The second step 
is to Pade-approximate the Taylor expansion of these 
functions. The third step is to use these Pade trans- 
forms in the definition of the energy E (Eq. ^ and Pade- 
approximated flux (Eg. [^0)1 . The last step is to construct 
either the Pade-approximated waveform as in Sec. llVl or 
the EOB waveform as in Sec. |Vl Schematically, the sug- 
gested procedure is summarized by the following map: 
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m, rm 
1 Jn 



£^(e"),f^(/r) 



(3) 



Our notation is to denote by T^{x) the Pade approx- 
imant of a fc-th order Taylor series Tk{x) with an m-th 
order polynomial in the numerator and an n-th order 
polynomial in the denominator such that m + n < k, i.e. 
the Pade approximant of ek{x) is e^{x). 

In Section ini we compare the 3PN Taylor series of the 
energy function to its possible Pade approximants using 
the intermediate energy function e(x), as suggested by 
Damour et al. [s^]- We compute the last stable orbit 
frequency, defined as the frequency for which the energy 
reaches a minimum as a function of frequency, and also 
the poles of the energy in the complex plane correspond- 
ing to each possible Pade approximant. The large vari- 
ation of last stable orbit frequency and poles does not 



PN order 


F„ 
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F„ 




0.0 


1.000000 


1.000000 


1.530011 


1.530011 


0.5 


1.000000 


1.000000 


1.000000 


1.000000 


1.0 


0.851547 


1.000000 


0.772866 


0.602534 


1.5 


0.952078 


0.911487 


1.005361 


0.887757 


2.0 


0.944193 


0.928720 


0.940013 


0.937227 


2.5 


0.931939 


0.936461 


0.925444 


0.938929 


3.0 


0.941025 


0.939366 


0. 945405 


0.939502 


3.5 


0.939726 


0.939399 


0.938991 


0. 938082 


4.0 


0.939208 


0.939363 


0. 939048 


0.939471 


4.5 


0.939745 


0.939719 


0.939979 


0.939516 


5.0 


0.939601 


0.939653 


0.939526 


0.939684 


5.5 


0.939605 


0.939623 


0.939616 


0.939621 



TABLE II: Convergence of the Taylor series and its Fade 
aprroximants of the flux in the test particle limit at v = 0.2 
{x = 0.04). The four flux functions Fn, FZ^', F„ and F^+' 
are given in Eqs. 1141 1151 1181 and [20] respectively. Even in the 
test mass limit, the Fade approximant of the flux fails to con- 
verge faster that its 5.5 PN Taylor series at a relatively small 
value of u = 0.2. After 12 terms, only about 4 or 5 signifi- 
cant digits seem reliable for the Taylor expansions and their 
Fade approximants. The lack of improvement in the conver- 
gence of the Fade approximants should be contrasted with 
the example in Table H) 



suggest good convergence of the Pade-approximated in- 
termediate energy function e{x). The energy function 
E{x) is strongly dependent on the choice of the Pade 
approximant of e{x). Accordingly, the Pade waveform 
will also be strongly dependent on the choice of the Pade 
approximant. 

In Sect ion Uni we present two possible methods for cal- 
culating the Pade approximant of the flux function. The 
first method simply takes the Pade approximant of the 
Taylor scries trea ting the logarithmic contribution as con- 
stant. Following [3J], the second method adds a pole to 
the Taylor series, factors out the logarithmic contribution 
to the series, and then computes the Pade approximant 
of the resulting Taylor series. We test the convergence 
of the Pade approximant for both methods versus their 
Taylor series. We find that the Pade approximants of the 
fiux do not converge any faster than their Taylor coun- 
terpart. 

A simple example that illustrates the problem is shown 
in Table [ill There we compare the partial sums of the 
Taylor series for the flux with the corresponding Pade 
approximants in the test mass limit. The four flux func- 
tions En, fZ^\ En and are given in Eqs.[T4l[15l[l8| 
and [201 respectively. Even for a relatively small value of x, 
namely x = 0.04 {v = 0.2), the Taylor series is converging 
very slowly. After 12 terms, only about 4 or 5 significant 
digits seem reliable. Moreover, the Pade resummation 
shows very similar behavior; there is no improvement in 
the convergence. We will return to this example in Fig.[3l 

In Section IIVI we generate all the possible Pade wave- 
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forms as suggested by Damour et al. [3J| corresponding 
to 3 and 3.5 PN order. The waveform approximation 
requires the choice of a pole. We use the only physi- 
cal pole, found from the 2PN Fade-approximated energy 
El. We also use the last stable orbit from the 3PN en- 
ergy Taylor series E3. The results are not very sensi- 
tive to this choice. We compare the Fade waveforms to 
a 15-orbit numerical waveform in the equal mass, non- 
spinning quasi-circular case [67j . The phase difference 
in these comparisons ranges between 0.05 and a few ra- 
dians for well-defined Fade approximants (not having a 
pole in the frequency domain of interest) when the match- 
ing of the numerical and Fade waveforms is done at the 
gravitational wave frequency AIui = 0.1 [63|. None of 
the Fade waveforms agrees with the numerical waveform 
better than the Taylor series T4-3.5/3.0FN, which has 
an error of 0.02 radians. (We identify post-Newtonian 
approximants with three pieces of information: the label 
introduced by [35] for how the orbital phase is evolved; 
the FN order to which the orbital phase is computed; 
and the FN order at which the amplitude of the wave- 
form is computed. See Ref. [g^ for more details.) Our 
conclusion is that the Fade approximant might be help- 
ful in suggesting fitting formulas but it does not provide 
a more rapidly convergent method. Note that the Fade 
transform also fails to accelerate the convergence of the 
T2, T3 and Taylor series (see Refs. [H, [g^ for the 
definition of these Taylor series). 

In Section |Vl based on the results of the previous sec- 
tions, we design a simple EOB model (closely related to 
the ET EOB model of Ref. ^) using the Taylor series of 
the flux. We add one unknown 4FN term that we fit for 
by maximizing the agreement between the EOB model 
waveform and the numerical waveform. The model does 
not require adding a pole to the flux, nor an a priori 
knowledge of the last stable orbit from the energy func- 
tion. This simple EOB model, with only one parameter 
to fit for, agrees with the numerical waveform to within 
0.002 radians (3 x 10^** cycles). (This is six times smaller 
than the claimed numerical accuracy of [39], smaller by 
an even larger factor than the claimed numerical accu- 
racy of [4^, and 25 times smaller than the gravitational 
wave phase uncertainty of the numerical waveform. See 
Table III in Ref [g^] for more details.) This model agrees 
with the numerical waveform better than any previously 
suggested Taylor, Fade or EOB waveform. 



where mi, m2 are the masses of the bodies. The total 
relativistic energy function E^tot is related to the post 
Newtonian energy function E{x) through 



Etot{x)^M[l + E{x)] 



(5) 



where M is the total mass (M — mi + ni2). Solvin g fo r 
E{x) in terms of e{x) using Eqs. ^ and ([5]), we get 34] 



E{x) = {1 + 21/ + <x) - 1] - 1 , (6) 

where the symmetric mass ratio is = mim2/M^. The 
orbital energy function E(x) is known as a Taylor series 
Ek up to 3FN order as a function of x and v [15| 
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Using the above equations, we compute the Taylor series 
expansion, efe(x), of e{x) up to 3FN order: 



e3PN(a;) 



1 -iK 

1 - (1 + -^.)a; - (3 - -i.)x2 



9-f — (-17236 -|-6157r2):/ 
288 



103 

1 

36 



81 



(8) 



In the test mass limit — + 0), the exact function e{x) 
coincides with the Fade approximant (x) of its Taylor 
expansion in Eq. ([5]) 



eix: V 



0) 



1 -4a; 
1 -3a; 



(9) 



This quantity has a pole at Xpoic = 1/3. The orbital 
energy is then 



E{x\v 0) = 



^l-x'-^-X 



1 -3.T 



(10) 



II. ENERGY FUNCTION 

Damour et al. [s^] introduced a new energy-type func- 
tion e(a;), where x is the FN frequency related parameter. 
This assumed more "basic" energy function eix) is con- 
structed out of the total relativistic energy £'tot(a;) of the 
binary system. Explicitly 



e{x) 



Etot 



2mim2 



- 1 



(4) 



1 - 6a; 



and it derivative is 

dE{x; _ 

^ "''2(l-3a;)3/2 

The last stable orbit occurs where 
ax 



(11) 



(12) 



so in the limit i> ^ the last stable orbit is at exactly 
^iso = 1/6. On the grounds of mathematical continu- 
ity between the test mass limit z/ — > and the finite 
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FIG. 2: Post Newtonian Energy at 3PN and its Fade approx- 
imants for the case v = \|^:. The plot includes the high value 
of x\^o = 0.36, the numerical data available is at x = 0.16. 
The plots of i?0) ^'i) E\ and E% vary significantly although 
they all correspond to the 3FN Taylor series of the energy 
function. E\ is very diff'erent from the other functions, which 
suggests a poorly convergent Fade approximant. 

mass ratio case, Damour et al. 34] argued that the ex- 
act function eix) should be meromorphically extendable 
in at least part of the complex plane and should have a 
simple pole on the real axis. They suggested that Pade 
approximants would be excellent tools for giving accurate 
representations of functions having such poles. 

Once we know the Taylor series of the new energy func- 
tion, efe(x), we compute its Pade approximant e™(a;), 
with m + n < k. The Pade-approximated energy, 
is obtained by replacing e(x) in Eq. ^ with e™(a;). In 
the equal mass case {v = 1/4), we can define several 
Pade approximants of ek{x). The most interesting Pade 
approximants have a maximal sum of their indices, since 
they should be closest to the unknown exact function if 
the Pade resummation is converging. In Fig. [21 we show 
a plot of the PN energy function E3P^{x) and its Pade 
approximants El, Ef, E\, E^ and Eq as a function of x. 

Although the Pade approximants of the energy are of 
maximal order, they differ significantly. Good conver- 
gence of the Pade approximants requires good agreement 
between approximants of the same order n + m, if there 
is no pole in the region of interest {0 < x < 0.4). For ex- 
ample, there is no a priori reason why one should prefer 
either Ef or E2 ■ Although both have the same order and 
are equally close to the diagonal, the difference between 
these functions is quite large. 

In Table Hill we compute the locations of the poles and 
the last stable orbits for all of these Pade approximants. 
The ill-convergence of the Pade transform is again seen by 
looking at the variation of the last stable orbit positions. 
In Table Hill for example, xiso of Ef differs by about 8% 



from xiso of i?2- Moreover, for finite v, the poles are all 
complex or not in the interval [0, 1] except for the case 
a^poic = 52/109, corresponding to the Pade-approximated 
energy El. There is no reason why this should be the 
"exact" pole that should be used in the formalism, since 
none of the third-order Pade approximants of the 3PN 
energy has a physical pole. 

In summary, using Pade approximants for the energy 
function in the equal mass case does not seem to pro- 
vide any benefit. The differences between the various 
Pade approximants of the energy are large. The quanti- 
ties Xpoio and Xiso do not show any regular behavior that 
could be a sign of a physical pole that could be found by 
using the Pade transform. 

III. FLUX FUNCTION 

The general form of the PN flux at order N is: 

F(v) = X , (13) 

5 

where the normalized flux is a Taylor expasion in v 
with logarithmic terms 

2N /2N \ 

fn{v) = J2 + J2 ' (14) 

fc=0 \fe=6 / 

where the post-Newtonian coefficients Ai and Bi are 
functions of the mass ratio parameter v. They are given 
in the test mass limit in Ref. [69ll and in the equal mass 
quasi-circular case in Ref. [l5|. The flux series has a 
logarithmic contribution starting at 3PN. Pade approx- 
imants, however, are well defined only for pure polyno- 
mials. Two possible methods are therefore used to com- 
pute the Pade approximant of the flux. The first method 
simply treats the logarithmic terms as constants and re- 
sums the series as a pure polynomial such that the Pade- 
approximated flux F„ is 

The second method, suggested by Ref. [ll], defines a 
new flux function / by adding a pole, factoring the log- 
arithmic terms from the series, and finally computing 
the Pade approximant of the pure polynomial. Since 
we would like to check the convergence of the Pade- 
approximated flux versus its Taylor series, we sketch the 
definitions of the various functions involved. According 
to Ref. |3J] , two ideas are needed for a good representa- 
tion of the analytic structure of the flux. First, since in 
the test mass limit F is thought to have a simple pole at 
the light ring [4^, one might expect it by continuity to 
have a pole in the comparable mass case. This motivates 
the introduction of the following factored flux function, 

/(«;^): 

f{v;>^)^(l ^)f(«;^^), (16) 
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Energy 


^polc 




-EsPN 




0.254 


El 


52/109 = 0.477 


0.199 


Eo 




0.262 


El 


-4.41 


0.261 


El 


0.170 ±0.757i 


0.285 


El 


0.044 ±0.50H,-0.696 


0.363 



TABLE III: Values of the poles and last stable orbit (Iso) of 
the energy for the case v — 1/4. The poles a;poie and last 
stable orbit frequency of the function E^{x) depend signif- 
icantly on which Pade approximant is constructed from the 
Taylor series Ckix). The only physical pole is a;poiG = 52/109, 
which is at a larger value than the pole in the test mass limit. 
The position of the last stable orbit also varies significantly. 



where Wpoio is the pole of the Pade-approximanted energy 
function used. 

Second, the logarithmic term that appears in the flux 
function needs to be factored out so we can use the stan- 
dard Pade transformation. After factoring the logarith- 
mic terms out, the flux function / becomes 



fn{v;v) = 



' 2N 



1 + log ■ 




(17) 



where the coefficients Ik and fk are given in Ref . [3J] , and 
v\so is the velocity of the last stable orbit of the Pade- 
approximated energy. Then the Taylor series of the flux 
with a pole is defined as 



Fn{v; v) 



^ fn{v\v) 

1 - v/Vpo\c{v) 



(18) 



The Pade approximant of the intermediate flux function 
f{v) is deflned as 



1 + log ■ 



Vis 



■27V 

E 



(19) 



where v\so{e^]v) denotes the last stable orbit velocity 
for the Pade approximant P™[e(x)]. Finally, the cor- 
responding Pade approximant of the flux F{v) is given 
by 



1 - v/vpoic{e'^]v) 



(20) 



where Wpoio(e™;i/) denotes the pole velocity defined by 



Flux for the test mass case 



The exact gravitational wave luminosity F is not 
known analytically in the test particle limit. It has 
been computed numerically by Poisson ^68^. The post- 
Newtonian expansion of the flux is known in the test mass 
limit to 5.5PN order [6^. This allows us to test the rate of 
convergence of the Taylor series of the normalized flux Fn 
fEg.fH)) and its Pade- approximant F„ constructed treat- 
ing the logarithmic term as a constant (Eg. [H)) . We also 
test the convergence of the flux function F„ (Eq. [T8|) and 
its Pade approximant P™ (Eq. EO)) . These convergence 
tests use the known values VpoXc = 1/ %/3 and v\so = 1/ V6 
for the test mass limit as discussed in Sec. [ill 



In Fig. [31 we test the convergence of the various flux 
functions at the velocity value v = 0.2. The four flux 
fimctions P„, P™'*''^, Fn and P™^*^ are given in Eqs. fT4l 
[TSlfTSl andl^resDectivelv. We use the Pade approximant 
along the diagonal P^^*^ where e = or 1. The rates 
of convergence of the Taylor expansion and its Pade ap- 
proximant are nearly equal for the two methods, whether 
or not we include a pole. As the PN order increases, the 
Taylor series and its Pade approximant alternate in which 
provides a better fit to the numerical data for the flux. 
For example, at 2PN order the Taylor flux with a pole 
(Eq. \TE\i flts the numerical data the best. At 2.5 and 3 
PN order the Pade approximant of the flux P™ (Eq. 
fits the numerical data the best, while at 3.5 and 5PN or- 
der the Taylor series of the flux (Eg. [H]) is the best. At 
5.5PN the Pade approximant of the flux (Eq. [20]) gives 
the best agreement. The results are similar for other 
values of v. No method has the best convergence rate. 



According to Pade theory, the convergence of the Pade 
approximant is best along the diagonal, but it is equally 
good along the off-diagonal terms if no pole exists in the 
region of interest (i.e. no zeroes appear in the denom- 
inator of the Pade approximant.) For this reason, we 
show the error between all the possible maximal Pade- 
approximated fluxes P„ (Eq. [T4|) and the numerical 
flux for three values of v (= 0.2, 0.25, 0.35) {x = 0.04, 
0.06, 0.12) in FigH The 5.5PN Taylor series, denoted by 

Pq , fits the exact numerical data better than the Pade 

approximants F^ ,Pg,Pg,Pg. In the other cases, the 

Pade approximants provide a better agreement (i.e. Pj^q, 

P3, P4 and P5) for the three values of v. This suggests 
that the Pade approximation should only be used to sug- 
gest a fitting formula for the numerical data, since there 
is no internal self-consistency in the agreement. The off- 
diagonal approximants do not show any regular pattern 
of convergence to the numerical data nor are they better 
than the Taylor series. 
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FIG. 3: Convergence of the flux approximations in the test 
mass hmit for v = 0.2. The four flux functions F„, F™^", F„ 
and are given in Eqs.[ll[l5l[l8]and[20]respectively. The 

Fade approximants do not converge faster than their Taylor 
series counterparts. The Fade and Taylor series alternate at 
providing the best agreement with the "exact" data as the 
FN order increases. Contrast the behavior here with Fig. [1] 



FIG. 4: Error between maximal Fade approximants of the flux 
F fEg. llSf) and the numerical flux in the test mass limit at n = 
0.2,0.25,0.35. The 5.5FN Taylor series, denoted by f", fits 



the exact numerical data better than the Fade approximants 
Fi , Fq, Fg, Fq. In the other cases, the Fade approximants 
provide a better agreement (i.e. Fig, F2, -F3, -F4 and F^, F^ 
and Fii). 



B. Flux for the equal mass case 

For binaries of comparable mass on a quasi-circular 
orbit, the flux is known only to 3.5PN order ^1^. In 
Ref. for a quasi-circular non-spinning binary, the 
numerical flux was computed by integrating the spin- 
weighted spherical harmonic components of the Weyl 
scalar "^4. The numerical flux data we use in this paper 
was provided by Harald P. Pfeiffer and Michael Boyle. 
The estimated error in measuring the flux data was about 
0.2%. The velocity range for the simulation was from 
V ~ 0.26 (x ~ 0.06) to w ~ 0.4 {x ~ 0.16). 

In the equal mass case, we cannot do an accurate con- 
vergence test early in the evolution as in Fig. [3] for two 
reasons. The flrst reason is the "junk radiation" (noise 
early in the evolution from imprecise initial data) during 
the first few orbits. The second reason is the inability 
to accurately define the numerical flux as a function of 
the orbital frequency of the binary. The numerical nor- 
malized flux is computed as a function of ^^22/2, where 
W22 is the wave frequency of the /122 mode. Instead, in 
Table IIVI we compare the convergence of the four flux 
fimctions F„, F™^^, F„ and (deflned in Egs. [Til 

[TSl [T8l and [20] respectively as a function of PN order) 
for V = 0.2 {x = 0.04), Wpoio = 0.69 (xpoic = 52/109) 
and viso — 0.50 (a;iso — 0.254). We use the last stable 
orbit frequency corresponding to the 3PN Taylor series 
of the energy and the pole corresponding to El. The 
convergence does not depend on these values although 
the flux values listed in Table IIVI do depend somewhat 



on the values of t^poie and viso- We choose a medium ve- 
locity {v = 0.2) to make the rate of convergence clear. 
At 3.5PN order, all four flux functions agree to 2 signif- 
icant figures. However after seven terms, Fn converged 
to 3 significant figures, F^^*^ converged to 4 significant 
flgures, while F„ and F^'^'^ converged to 2 significant 

flgures. The flux function F^ converged to one addi- 
tional signiflcant flgure over Fn, however F^ cannot 
reliably be considered more accurate than Fn because it 
converges to a slightly different value. The Pade approx- 
imants do not seem to converge to a larger number of 
signiflcant flgures than the Taylor flux function Fn- 

In Fig. [51 we plot the numerical normalized flux Fnr, 
the 3.5PN flux ^"3.5 and the maximal Pade-approximated 
flux functions F|, F^, F^, F^ and F^ {= F7). Although 
^"3.5 diverges from the numerical flux early at t; ~ 0.26, 
it still fits the numerical data better than F^, Ff and Fq. 
The quantity has a pole and fails to capture the nu- 
merical flux behavior completely. The quantity Fq is by 
deflnition the Taylor flux with a pole, F7. This function 
shows that adding a pole to the Taylor expansion of the 
flux i^3.5 degrades the flt with the numerical flux. More- 
over, the numerical flux does not suggest the existence of 
a pole at a large velocity {v ~ 0.69); it starts to decrease 
to at u ~ 0.4. Adding a pole does not seem a useful 
idea in this case at least. On the other hand, i^l ^^'^ ^4 
are a better flt to the numerical data during most of the 
velocity range of the 15-orbit data. The flux function _F| 
is especially a good flt to the numerical normalized flux 
at high velocities. However, even though and F^ are 
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F„ 


' m 




^ m 


0.0 


1.000000 


1.000000 


1.407582 


1.407582 


0.5 


1.000000 


1.000000 


1.000000 


1.000000 


1.0 


0.822381 


1.000000 


0.749987 


0.353292 


1.5 


0.922912 


0.886577 


0.963887 


0.865262 


2.0 


0.922745 


0.905792 


0.922678 


0.910047 


2.5 


0.904387 


0.910595 


0.896904 


0.912033 


3.0 


0.913204 


0.912261 


0.916323 


0.912613 


3.5 


0.913314 


0.912223 


0.913275 


0.911492 



TABLE IV: Flux convergence in the equal mass case for 
V = 0.2 {x = 0.04), Upoic = 0.69 (a;poic = 52/109) and 
Visa = 0.50 (a;iso = 0.254). The four flux functions Fn, -F'm^'^, 
Fn and are given in Eqs. [Til [TSl [TSl and [20l respectively. 

At 3.5PN order, all four flux functions agree to 2 signiflcant 
figures. After seven terms, F„ converges to 3 significant fig- 
ures, -Fm^' converges to 4 significant figures, while Fm and 
F^^'^ converge to 2 significant figures. 




a good fit to the numerical flux during the last 15-orbit 
inspiral before merger, there is no guarantee that this is 
true at low velocities. 



FIG. 5: Normalized fiux for an equal mass nonspinning bi- 
nary. We plot the numerical fiux Fmr, the 3.5PN flux fs.s 
and the maximal Pade-approximated flux functions Ff, F^, 
Fi, F^ and (= F7). The early noise is caused by the junk 
radiation. 



IV. FADE WAVEFORMS 

The construction of the post-Newtonian waveforms re- 
quires solving the post-Newtonian equations describing 
the motion of the binary and the generation of gravita- 
tional waves. Substituting the orbital evolution predicted 
by the equations of motion into the expressions for the 
waveform would not generate waveforms accurate enough 
for matched filtering in detecting gravitational waves [70| ■ 
To compute the waveform at 3PN order, it is necessary to 
solve the equations of motion at 5.5PN order, because the 
radiation reaction contributes to the equations of motion 
starting at 2.5PN order. However, for a non-spinning bi- 
nary of equal mass and on a circular orbit, accurate wave- 
forms at 3PN order can be constructed under two further 
assumptions. The first assumption is that the binary fol- 
lows a slow adiabatic inspiral. The second assumption is 
that of energy balance between the orbital binding en- 
ergy and the energy emitted by the gravitational waves, 
where the energy balance equation is defined as 



dt 



(21) 



The procedure of constructing the standard Pade 
waveforms [s^l is similar to one used to construct the 
TaylorTl waveforms in Ref. [33,[63|. The main difference 
is the use of Pade approximants of the energy and flux to 
compute the orbital phase, as described in Sees. |TT] and 
mil instead of their Taylor expansions. The orbital phase 
used in the Pade waveforms is obtained by numerically 
integrating 



'dt 



The fraction on the right side of Eq. ([22)1 is retained as 
a ratio of the Pade approximants of the post-Newtonian 
expansions, and is not expanded further before numerical 
integration. The waveform is produced by substituting 
the orbital phase into the spherical harmonic mode /i22 
of the post-Newtonian waveform, which is known up to 
3PN order [H, [H, Is^, HH . 

Given the expressions for the Pade-approximated en- 
ergy and flux in Sections HIl and IIIIl and the Taylor series 



of the waveform amplitude 



there is still 



dEfldVL 



(22) 



a set of choices that must be made in order to produce 
a Pade-approximated waveform that can be compared to 
our numerical waveform. These include 

1 . The Pade approximant of the orbital energy, Ef . 

2. The flux function and its Pade approximant F^. 

3. The velocity of the pole and the last stable orbit, 
Wpoic and v\so- 

4. The PN order through which terms in the waveform 
amplitude are kept. 



A. Procedure 

We consider numerical gravitational waves extracted 
with the Newman-Penrose scalar ^4, using the same pro- 
cedure as in [7l|. To minimize gauge effects, we compare 
its (2, 2) component extrapolated to infinite extraction 
radius according to Ref. [67| . The extracted waveform is 
split into real phase (j) and real amplitude A, defined by 
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FIG. 6: Phase difference between the 3 and 3.5 PN Pade 
approximated and numerical waveforms matched at the wave 
frequency Mu = 0.1. We use the Pade-approximated flux 
FY!Lm (Eq. [20]) and energy Ef . We include in the figure the 
waveforms using the Pade-approximated flux using m — 
— 1. There is no entry for m = A since the Pade-approximated 
flux has a pole in the frequency range of the simulation. 



Ref. [63] as 

^-f (r,t) = A(r,Oe"'^("'*^ (23) 
The gravitational-wave frequency is given by 

The spherical harmonic component (2,2) of is then 
compared to the numerically twice-differentiated post- 
Newtonian expression of /122 , ^22 , as in Ref. [63| . Follow- 
ing [73, [tJ], the matching procedure needed to set 
the arbitrary time offset tg and the arbitrary phase offset 
00 is done by demanding that the PN and NR gravita- 
tional wave phase and gravitational wave frequency agree 
at some fiducial frequency ujm- 

B. Results 

In this section, we compare the numerical waveform 
to the Pade waveforms corresponding to the 3.5 PN or- 
der of energy and flux using the 3PN Taylor series of 
the post-Newtonian amplitude A22 ■ The energy and flux 
functions used are those suggested by Ref. f34!|. We do 
not generate all possible waveforms using different Pade 
approximants of the energy or the flux at low PN orders, 
since all these resummed series showed no improvement 
in the convergence rate. 

As introduced in Sec. [TTl we use the Pade- 
approximated energy Eg, Ef, E\ and corresponding 



to the PN Taylor series of the energy, and the Pade- 
approximated energy, E\ , corresponding to its 2PN Tay- 
lor expansion. For the flux, the diagonal Pade approxi- 
mant F3 is used in addition to all possible Pade approx- 
imants of flux at 3.5PN order, F^tm: where < to < 7 
as described in Sec. IIIII 

The Pade-approximated flux has two parameters, v\so 
and Wpoio as discussed in Sec. IIIII The value fpoic — 
52/109 is used. We also tested varying the pole loca- 
tion, but found that we could not improve the agreement 
significantly. 

From Table IIIII any value of the velocity of the last 
stable orbit could be used. We use the 3PN value v\so = 
0.254 and also use v\so = 0.199. The latter is used when 
the Pade approximant E\ is employed in the construction 
of the waveform. In the remaining cases, we use v\so — 
0.254 since it is quite close to the estimates from other 
Pade approximants of the energy. The effect of changing 
the value of v\so is not significant compared to changing 
the order of the Pade approximant for the energy or the 
flux. 

To do the comparison, we match the Pade- 
approximated and numerical waveforms at the wave fre- 
quency Mlo = 0.1. Then we measure the maximum 
phase difference between the numerical waveform and 
each of these Pade waveforms during the inspiral when 
the numerical wave frequency is between Mlo = 0.035 
and Mid = 0.1 (as in the upper panel of Fig. [7]). Our 
results are summarized in Fig. [SI which shows the phase 
differences for each of the Pade approximants of energy 
Ef and flux F^^. On the same Figure, we include phase 
differences for the waveforms generated using the Pade- 
approximated flux under the to = — 1 entry. 

When El is used, the phase error ranges between 2 
and 5 radians as to increases from —1 to 7. Using all the 
possible Pade approximants of the 3PN energy, the esti- 
mated phase difference ranges from 0.05 to 2.5 radians. 
Using the Taylor series with a pole (m = 7) resulted in a 
large phase difference ranging between 1 and 1.5 radians. 
The diagonal Pade term F^ of the flux generates similar 
phase differences, ranging from 0.06 to 0.2 radians as the 
Pade order of the energy changes. 

The Pade-approximated waveforms do not fit the nu- 
merical data better than the waveforms using the Tay- 
lor expansion of the flux. Although the Pade-waveforms 
along the diagonal have a phase difference less than 0.25 
radians, none of these waveforms fits the numerical wave- 
forms better than TaylorT4 at 3.5PN order as shown in 
Ref. [131 . Moreover, the dependence of the phase differ- 
ence on the Pade order suggests that there is no reason 
why it should help in estimating the parameters better in 
data analysis. This is as expected from the poor conver- 
gence of the the Pade approximant of the ffux discussed 
in Sec. [nH 

The Pade resummation techniques were also tested on 
the Taylor series for the amplitude, and they showed no 
improvement in the convergence of the series. In adi- 
tion, none of the tests that were performed on the Pade 
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resummed Taylor series of the T2 and T3 waveforms 
showed a faster convergence rate. In fact, there is no 
improvement in convergence for any Taylor series in the 
PN approximation that we have investigated. 



V. SIMPLE EOB MODEL 

We have described the failure of the Pade resumma- 
tion techniques to accelerate the convergence of any PN 
Taylor series, the absence of any signature of a pole in 
the flux in the equal mass case, and the erratic pattern 
of agreement between the Pade waveforms and the nu- 
merical waveform. It seems one might as well simply use 
the Taylor series at all steps of computing waveforms. 
Also it does not seem that the parameters Vpoie and f iso 
are useful. In this section, we show how to get good 
agreement with the numerical waveform by using a sim- 
ple EOB model. The only parameter we introduce and 
fit for is an unknown 4PN contribution to the flux. 



A. EOB waveforms 

The EOB formalism (53| is a non-perturbative ana- 
lytic approach that handles the relative dynamics of two 
relativistic bodies. This approach of resumming the PN 
theory is expected to extend the validity of the PN results 
into the strong- field limit. The procedure for generating 
an EOB waveform follows closely the steps in Sec. lIVI In- 
stead of using the energy balance equation, we compute 
the orbital phase by numerically integrating Hamilton's 
equations. The EOB waveform is generated by substitut- 
ing the orbital phase into the waveform amplitude A22 at 
3PN order. The two fundamental ingredients that allow 
computing the orbital phase are the real Hamiltonian H 
and the the radiation reaction J-'a.. 



B. Hamilton's equations 

In terms of the canonical position variables r and (j) and 
their conjugate canonical momenta Pr and p^, where r is 
the relative separation and 4> is the orbital phase, the real 
dynamical Hamiltonian is defined as [53 |: 



1 



H ^ -^l + 2:y {Heob -1). 

V 



(25) 



where 



ff^OB = J a(i + ^ + § + 2.(4-3.)^), (26) 



The Taylor series of the A function is replaced by its Pade 
approximant ^3 . Here the Pade approximant is not used 
to accelerate the convergence of the Taylor expansion of 
A. Instead, it leads to the existence of a last stable orbit 
(see Ref. [39[ and references therein). Otherwise, the 
EOB Hamiltonian is non-physical for the last few orbits; 
the orbital frequency stays nearly constant for several 
orbits before merger. For the B function, the Taylor 
expansion suffices: 



B = 



1 - 



2(3i/- 26)- 



(28) 



Then Hamilton's equations of motion are given in the 
quasi-circular case by 



dtr 

dtPr 



(29) 
(30) 

(31) 
(32) 



where Tij, is the radiation reaction in the cj) direction rep- 
resenting the nonconservative part of the dynamics. In 
Eq. [32I dij,H = since H is independent of 4>. The radi- 
ation reaction is deduced from the post-Newtonian flux 
as in Refs. Eil,!!!,!!! 



J^d, — 



F + F.v^ 



(33) 



In this equation, we have introduced an unknown 4PN 
flux term, the only parameter that we fit for in this 
EOB model. 



C. Initial conditions 

To integrate Hamilton's equations, we need appropri- 
ate initial conditions for a quasi-circular orbit. Refs. |4ll . 
m, HI] indicate how to define some "post-adiabatic" ini- 
tial conditions. However, these initial conditions do not 
generate an orbit with as low an eccentricity as the nu- 
merical simulation, roughly 5 x 10^^. At a given radius r, 
starting from the post-adiabatic initial conditions of Pr 
and we therefore reduce the eccentricity iteratively 
in two steps. The first step includes evolving Hamilton's 
equations in the conservative regime {J^ = 0) and it- 
eratively changing the value of Prj, until the eccentricity 
measured from the evolution of the orbital separation is 
of the order 10~^. The second step is based on evolving 
the nonconservative Hamilton's equations with the 4PN 
flux and iteratively changing the Pr momentum until the 
eccentricity is again of the order 10~^. This circulariza- 
tion procedure is repeated as we iterate to maximize 
the agreement between the waveforms. 



and where the radial potential A function is defined as 
the series 



A = 1 



2v 



/94 41 2\ 
ly~32^ Jr^ 



(27) 



D. Best Fit of Fs 

To find the best fit for Fg, we iteratively solve for the 
minimum in the phase difference between the numerical 
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(t-r*)/m 

FIG. 7: Phase and amplitude differences between the EOB 
waveform and the numerical waveform. After fitting for the 
best value of Fs, the phase difference is less than 0.002. The 
early noise is due to junk radiation at the early stage of the 
numerical simulation, r* is the tortoise coordinate defined 
in M. 



and EOB waveforms. The waveforms are matched as in 
Sec. II VI at the wave frequency muj = 0.1, and the phase 
difference is defined as the maximal phase difference dur- 
ing the inspiral phase up to the wave frequency muj = 0.1. 
We find a best fit value Fg = —333.75 corresponding 
to the initial conditions r — 17, = 0, Pr = —0.0008, 
P0 = 4.53235. A change of 1% in Fg changes the max- 
imal phase difference from less than 0.002 radians to 
about 0.01 radians. Note that without adding the fitting 
parameter Fg, the phase difference is about 1.7 radians 
during the 15-orbit inspiral. 



E. Results 

In the upper panel of Fig. [71 we plot the phase dif- 
ference between the numerical waveform and the EOB 
waveform computed using the 3PN Taylor series of the 
amplitude A22. The phase difference is less than 0.002 ra- 
dians after maximizing the agreement between the wave- 
forms in the region where muj < 0.1. The early noise is 
due to junk radiation at the early stage of the numeri- 
cal simulation as described in Sec II C of Ref. [g^l • The 
phase uncertainty in the simulation was estimated to be 
0.05 radians; See Table III in jl^. 

In the lower panel of Fig. [71 we plot the relative differ- 
ence between the amplitude of the numerical waveform 
and the EOB waveform. The EOB waveform amplitude 
does not fit the numerical waveform amplitude as well 
as the wave phase does. This is expected because the 
waveform amplitude is known to 3PN order only, and 
no free parameter in the amplitude was fitted for. The 



agreement between the amplitude of this EOB model and 
the numerical waveform is similar to the agreement be- 
tween the amplitude of TaylorT4 3.5/3.0 and the numer- 
ical waveform in Fig. 21 in |67|. 

This EOB model is a modification of the ET EOB 
model of Ref. [ill . It fits the numerical phase very well 
without using the Fade resummation techniques nor a 
pole in the flux. 

Even though we have found very good agreement be- 
tween the waveforms, these results only suggest that the 
EOB model is a very good fitting model. Moreover, hav- 
ing fit a particular waveform, there is no guarantee the 
model will have predictive power for a more general case. 



VI. CONCLUSIONS 

Convergence tests show that none of the Taylor series 
in the FN approximation, such as the energy or the flux, 
could be replaced by a Fade approximant that converges 
faster. Other attempts we tried to accelerate the conver- 
gence of these series also failed, as for example using the 
Levin method to accelerate convergence js^l- As a result, 
more reliable waveforms could not be constructed using 
a Fade resummation scheme. Moreover, the Fade wave- 
forms also do not fit numerical simulation data better 
than the Taylor waveforms. Thus, they do not seem to be 
better than the Taylor waveforms in building templates 
for waveforms. This conclusion is independent of the 
Fade approximants used to test the convergence. Taking 
for example the sub-diagonal Fade approximant does not 
show any improvement in the convergence rate. In addi- 
tion, this conclusion is independent of the numerical data 
we used. We can simply take the highest FN order of the 
Taylor series or the Fade approximant and use it as the 
"exact" value of the function to test the convergence at 
low frequency. 

Based on the dependence of the flux on the velocity 
in the equal mass case, we do not find it helpful to add 
a pole to the flux. Therefore, we recommend using Tay- 
lor series instead of the Fade approximant to generate 
waveforms both in the time and frequency domains. The 
simple EOB model used in this paper agrees with the 
numerical data very well; the phase difference during the 
inspiral is much less then the estimated phase uncertainty 
in the numerical data. This model does not use Fade ap- 
proximants or poles except in one place to enforce a last 
stable orbit. Since Fade approximation does not acceler- 
ate the convergence of any FN Taylor series, there is no 
reason why it should estimate parameters better in data 
analysis of waveforms. 
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